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Abstract 

In this paper we propose a solution to the problem of moment matching with preservation of the port Hamiltonian 
structure, in the framework of time-domain moment matching. We characterize several families of parameterized port 
Hamiltonian models that match the moments of a given port Hamiltonian system, at a set of finite interpolation 
points. We also discuss the problem of Markov parameters matching for linear systems as a moment matching problem 
for descriptor representations associated to the given system, at zero interpolation points. Solving this problem yields 
families of parameterized reduced order models that achieve Markov parameter matching. Finally, we apply these results 
to the port Hamiltonian case, resulting in families of parameterized reduced order port Hamiltonian approximations. 
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1 Introduction 

IjTj'nrt Hamiltonian systems represent an important class of systems used in modelling, analysis and control of physical 
(-h systems, see e.g. |46M38j . These representations are used in lumped parameter system analysis and control stemming from 

4-3 e.g., mechanical systems, electrical systems, electromechanical systems, power systems. One of the main features of such 

2 systems is the passivity property, i.e., the internal energy of the system at some time is lower than or equal to the sum 

G of the stored energy in the past and the externally supplied energy in the interval between the past and the present (see, 

e.g., [46). Passivity is essential for stability analysis and controller synthesis represented by techniques such as Passivity 

^H Based Control (PBC) and Interconnection and Damping Assignment (IDA), see e.g., [37]. For instance, in the case of large 
£> power systems, consisting of synchronous machines interconnected through transmission lines, the passivity (dissipativity) 
^ analysis of its port Hamiltonian modelling is important for the stability assessment of such systems and for further control 

Xj design, see e.g., [33J [Tg]. However, physical modelling often leads to (port Hamiltonian) systems of high dimension, usually 

l/-v difficult to analyse and simulate and unsuitable for control design. Hence, model reduction is called for. In the systems 
•and control literature there are many results on model order reduction with preservation of properties and/or structure. 

^■LFor example, passivity preserving model reduction is discussed in e.g., [3j [45l [23 [HI CE3 and structure preservation is 

£C) considered in e.g., [T3ir25] for symmetric systems, [33J for mechanical systems, [HI [Ml 132] f° r P or t Hamiltonian systems. 

. .In the problem of model reduction moment matching techniques represent an efficient tool, see e.g. [HH 151 [31 |4"T1 [T31 I3T1 15] 

^* for a complete overview for linear systems. Using a numerical approach based on Krylov projection methods the (reduced 

k> order) model is obtained by efficiently constructing a lower degree rational function that approximates a given transfer 

; 1 function (assumed rational) . The low degree rational function matches the given transfer function at various points in the 

tt complex plane. If the interpolation points are at zero, then the Pade approximation problem is solved. If the interpolation 

points are finite, then the general rational interpolation problem is solved. In the case of multiple-input multiple output 

(MIMO) systems the problem is called tangential interpolation, i.e., finding an approximation that interpolates a transfer 

matrix at selected points along selected directions, see [15j \4\ for further details. If the interpolation points are at infinity, 

the problem is called partial realization and has been studied in e.g., |35] and references therein. 

Alternatively, for single-input single-output (SISO) systems, a system theoretic, time-domain approach to moment match- 
ing has been taken in [BJ. In short, the notion of moment of a linear, minimal system has been related to the unique 
solution of a Sylvester equation, see also, e.g., [T71 [IB], for previous results. Furthermore, the moments are in one-to-one 
relation with the steady-state response (provided it exists) of the given system driven by a signal generator, which "con- 
tains" the interpolation points. The moments have also been connected to the solution of a dual Sylvester equation, and 
shown to be in one-to-one relation with the well-defined steady-state response of the given system driving a (generalized) 
signal generator. Used for model reduction, the time-domain approach yields simple and direct characterizations of all 
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parameterized, reduced order models that match a prescribed set of moments of a given system at a set of finite inter- 
polation points. The classes of reduced order models that achieve moment matching contain subclasses of models that 
meet additional constraints, i.e., the free parameters are useful for enforcing properties such as, e.g., passivity, stability or 
relative degree, irrespective of the choice of interpolation points. 

Recently, in |49] the rational interpolation problem for linear port Hamiltonian systems has been addressed using Krylov 
projection methods, yielding reduced order models that match the moments of the given port Hamiltonian system at a 
set of prescribed finite or infinite interpolation points. Improved procedures for MIMO systems have been developed in 
[421 123) , where a near-optimal port Hamiltonian approximation that satisfies a set of tangential interpolation conditions is 
proposed. Furthermore, in 40, H2] the partial realization problem for port Hamiltonian systems has been considered. The 
procedure therein involves finding a change of coordinates such that the Hamiltonian becomes the square of the norm of 
the state vector, since a direct application of the Krylov methods does not yield a port Hamiltonian approximant. Then 
the Krylov projections are computed and applied to the system in the new coordinates, resulting in a port Hamiltonian 
model that matches the Markov parameters of the given system. Note that the Hamiltonian of the approximant is again 
the square of the norm of the reduced order state vector. 

In this paper we study the problem of computing low order approximations that match a set of prescribed moments at 
a set of finite or infinite points, of a given SISO, port Hamiltonian, linear system and preserve the port Hamiltonian 
structure, in the framework of time-domain moment matching, as described in the following problem. 

Problem 1 (General Approximation Problem). Given a linear, port Hamiltonian, SISO system, 

1. compute the families of parameterized reduced order models that match the moments of a given port Hamiltonian 
system at a set of finite or infinite interpolation points (in this case the models match a set of prescribed Markov 
parameters) . 

2. characterize the families of reduced order models which satisfy the following properties: the models match the 
moments of the given port Hamiltonian system and preserve the port Hamiltonian structure. In other words, from the 
classes of models that achieve moment matching we find the reduced order models that inherit the port Hamiltonian 
form. □ 

In the case of matching at finite interpolation points we obtain families of parameterized state-space, reduced order port 
Hamiltonian models that approximate the given port Hamiltonian system. All the reduced order state-space models 
share the same transfer function. In the SISO case, the state space parameters are used to enforce additional structure 
constraints such as diagonal Hamiltonian function, diagonal dissipation, etc. However, at the moment we are not able to 
determine which of the original variables and their meaning is retained in the reduced order model. In the MIMO case 
the free parameters can be used to define appropriate directions such that the reduced order models satisfy prescribed 
sets of tangential interpolation conditions. Regarding the computational aspect, we establish connections between the 
models from the families of parameterized port Hamiltonian reduced order models that achieve moment matching and the 
counterpart approximations obtained using Krylov projections, i.e., there is a relation between the free parameter and the 
Krylov projector. 

Since, to the best of the authors' knowledge, the existing moment matching techniques do not yield an a priori approxima- 
tion error bound and assuming the interpolation points are free parameters, we can use existing results from the literature 
to choose these parameters in order to improve the accuracy of the approximation, e.g., in the sense that the Hi norm of 
error decreases, although it is not the scope of the paper. 

We study the problem of Markov parameters matching, i.e., the partial realization problem, by extending the time-domain 
moment matching results to the case of interpolation points at infinity. We define the notion of moment for a class of 
linear, descriptor representations, associated to the transfer function of the given system, in terms of the (unique) solutions 
of generalized Sylvester equations and their dual counterparts. In particular, the Markov parameters of a given system 
are the moments of the associated descriptor realization at zero. Furthermore, we relate the moments to the steady- 
state response, provided it exists, of the descriptor realization driven by/driving signal generators. Performing model 
reduction we obtain several families of parameterized, (descriptor) reduced order models that match a set of prescribed 
moments of the descriptor realization associated to the transfer function of a given linear system. In particular, matching 
at zero yields classes of reduced order models that match the Markov parameters of the given linear system. As in the 
rational interpolation case, we establish relations between the parameter which define the family of the port Hamiltonian 
approximation that matches a set of Markov parameters and a Krylov projector, giving insight into the computational 
issue of the proposed solution. 

Finally, we apply these results to linear port Hamiltonian systems, resulting in families parameterized of state-space, 
reduced order port Hamiltonian models that match the Markov parameters of the given port Hamiltonian system. Note 
that the free parameters can be used to enforce additional physical structure to the approximant. We mention that, to our 
knowledge, there is no structured procedure to determine the number of interpolation points needed. However, based on 
[50] the number of interpolation points can be related to the order of the dynamics associated to the higher order Hankel 
singular values of the given system. To conclude, we mention that the scope of the paper is not to address computational 
issues, but to propose a system-theoretical based framework for port Hamiltonian reduced order modelling, consistent 
with the existing theory and suitable for future nonlinear extension. 



To give the complete picture, we compute several (equivalent) families of reduced order port Hamiltonian models, starting 
from the unique solution of a Sylvester equation and its dual counterpart. The families exhibit the same type of properties, 
without additional advantages or disadvantages, obtained however through different parameterizations. This argument is 
consistent with the rational Krylov projection modelling where the left and the right projection are dual to each other. 

The paper is organized as follows. In Section[5]we give a brief overview of the definition of moments and moment matching 
for linear port Hamiltonian systems, as well as of the family of parameterized reduced order models that achieve moment 
matching at a set of finite interpolation points. We also recall the procedures to obtain a port Hamiltonian approximation 
using Krylov projections for both SISO and MIMO systems. We present existing choices of interpolation points which yield 
accurate (in some sense) approximations. The section is completed with the formulation of the model reduction problem 
to be solved, a particular case of the [TJ for matching at finite interpolation points. In Section [3] we study the general 
problem of time-domain moment matching for a class of descriptor representations associated to the transfer function of 
a given linear system. In particular we show that matching the moments of the descriptor realization at zero is equivalent 
to matching the Markov parameters of the given system. Moreover, we obtain the classes of reduced order models that 
match a set of the Markov parameters of the given system. These models are then related to their Krylov projection 
counterparts. Based on Section [3J in Section U we discuss the problem of moment matching at a set of finite interpolation 
points, with preservation of the port Hamiltonian structure, and characterize the port Hamiltonian reduced order models. 
Furthermore, we give a necessary and sufficient condition for a reduced order model that achieves moment matching to be a 
port Hamiltonian model. We describe the families of parameterized state-space port Hamiltonian approximations, the free 
parameters of which can be used to enforce further properties. We also present a procedure that allows the computation 
of such a family of models. Based on the results from Section [31 in Section 14.21 we obtain the classes of reduced order, 
port Hamiltonian models that match the Markov parameters of a given port Hamiltonian system, proposing a procedure 
that allows the computation of such families of models. In Section [5J we give an example that illustrates the results. The 
paper is completed by a Conclusions section. 

Preliminary results have been presented in [25] . This paper constitutes a complete version of the aforementioned conference 
paper, with results presented, discussed and proven in detail. Furthermore, in this paper we develop new results based 
on the solution of the dual Sylvester equation resulting in an extended system-theoretic characterization of the moment 
matching for port Hamiltonian systems. In addition, a discussion of the multiple-input multiple-output case is given. 
Moreover, connections with previous work and results from the literature, especially Krylov projection-based modelling, 
are established and examples are given to illustrate the theory. 

Finally, note that this paper is a preliminary step to develop a time-domain moment matching based model reduction 
theory for nonlinear port Hamiltonian systems, see e.g., [26] for the first results in that direction. 

Notation. R is the set of real numbers and C is the set of complex numbers. C° is the set of complex numbers with 
zero real part and C~ denotes the set of complex numbers with negative real part. A* G c™ xm denotes the transpose and 
complex conjugate of the matrix A G C mxn . If A is a real matrix, then A* is the transpose of A. a(A) denotes the set 
of eigenvalues of the matrix A and denotes the empty set. nxi , G C nxi/ is the matrix with all elements equal to 0. If 
n = 1, then 0i x „ = [0,...,0] G U lxi/ and 0^ = 0* x „. S(t) denotes the Dirac (^-function. The triple (A,B,C) denotes the 
linear, time-invariant system described by the equations x = Ax + Bu, y = Cx, with input u(t) from a well defined set of 
inputs, output y{t) and state x(t). 

2 Preliminaries 

Let J G R nx ™ be a skew symmetric matrix and R G R" x ™, Q G R™ xn be two symmetric matrices. Consider the 
single-input, single-output, minimal, port Hamiltonian system 



x = (J- R)Qx + Bu, 
V = B*Qx, 



(1) 



where x(t) £ R™, u(t) G R, y(t) G R and B G R™. The Hamiltonian is H(x) = \x*Qx. The transfer function of system 
(JJ is given by K : C -> C, K(s) = B*Q(sI - (J - R)Q)- 1 B. Let s,eC,! = l, ..., v be such that s, <£ a{{J - R)Q). The 

moments of (JTJ at Si are Tio(si),T)x(Si), . . . , with rjk(si) = ^W d J s . Note that the moment rjk(si) represents the 

s= Si 

fc-th coefficient of the Taylor expansion of K(s) at Sj. Without loss of generality, throughout the paper we assume that 
fc = 0. 

2.1 Time-domain moment matching 

In this section we give a brief overview of a notion of moment in the time-domain setting, see [5] for a more detailed 
analysis. Based on this notion, families of parameterized reduced order models are developed. 



Consider the linear system (JTJ and let the matrices S £ W x " and L £ R lx ", and Q £ W x " and K £ W be such that the 
pair (L,S) is observable and the pair (Q,1Z) is controllable, respectively. Consider the Sylvester equation 

( J - R)QU + BL = ITS', (2) 

in the unknown II £ C nxv and its dual 

Qr = T(J~R)Q + KB*Q, (3) 

in the unknown T £ C vxn . Assume that a (A) n c(S) = 0. Since the system (JTJ) is minimal, the Sylvester equation © 
has a unique solution II and rank II = v. Similarly, if o~(A) D o~(Q) = 0, then equation ([3]) has a unique solution T and 
rank T = v. (see e.g., [TT]L 
Definition 1. 

1. Let (f) = [0i 02 •■• </>„] e <C lxv be such that 

= B*QIi. 

We call the moments of system (TTTjl at (7(5) the elements </>i, i = 1, ..., i/. The interpolation points are the eigenvalues 
of S, i.e., {81,82, ■■■, s„} = a(S). 

2. Let yj = [ipi <f2 ... ip v \* £ C be such that 

tp = TB. 

We call the moments of system (jTT)l at <t(Q) the elements </?j, i = 1, ..., z/. The interpolation points are the eigenvalues 
of Q, i.e., {si,S2,...,s„} = o-(Q). □ 

The following result establishes the connection between the notions of moment from Definition [TJ and the well-defined, 
steady-state response of the interconnection of {XJ with a signal generator. 
Theorem 1. [5| 0/ The following statements hold. 

1. Consider system (Qp and /ei £ 6e swc/i i/iai o~(S) C C . Assume that cr((J — R)Q) C C _ . Lei 

cj = 5w, w(0) ^ 0, (4) 

wra'i/i w(t) £ R." . Consider the interconnection of systems {Ip and Op, to£/i u = Lw, where L is such that the pair 
(L, S) is observable. Then the moments <f> of system {7P at o~(S) are in one-to-one relation with the well-defined 
steady-state response of the output y(t) = B*Qx(ty\of such interconnected system. 

2. Consider the system (QP and let Q be such that o~{Q) C C . Consider the system 

u = Qu + TZv, , , 

d = uj + Tx, [ ' 

where oj(t) £ W , d(t) £ W , v(t) £ R and 1Z is such that the pair (Q, 1Z) is controllable. Consider the interconnection 
between the system {7p and the system {5|], with v = y. Assume that <r((J — R)Q) C C - , a;(0) = 0, w(0) = and 
u{t) = 5{t). Then the moments tp of system fip at o~(Q) are in one-to-one relation with the well-defined steady-state 
response of the output d(t) . □ 

Based on Definition [TJ we define families of parameterized models of order v that match the moments of JJJ at the 
interpolation points {si, ..., s u } — o~(S) and families of parameterized models of order v that match the moments of ([TJ at 
{si,...,s„} = a(Q). 
Theorem 2. [JJ ^ The following statements hold. 

1. Let the pair (L, S) be observable and assume o~(A) n o~(S) = 0. Let £(t) £ W and consider the family of linear models 

i = (S-GL)S + Gu, () 

4> = B*Qnt;, y ' 

parameterized in G £ C v , where H is the unique solution of 0). Assume a(S — GL)Da(S) = 0. Let </> £ C lxi/ be the 
moments of @) at o~(S). Then (OJ) describes a family of reduced order models of f 11\ ), parameterized in G, achieving 
moment matching at o~(S), i.e., <f> = <f>. 

2. Let the pair (Q,1Z) be controllable and assume a(A) n c(<2) = 0- Let £(£) £ M? and consider the family of linear 
models 

y . f t=(Q-KH)Z + rBu, 

parameterized in H £ M lxi ' , where T is the unique solution of 0). Assume a{Q — 1ZH) n o~(Q) = 0. Let (p £ C lxu 
be the moments of Hj) at cr(Q). Then H\) describes a family of reduced order models of JJTp , parameterized in H, 
achieving moment matching at o~{Q), i.e., ip = (p. D 



1 By one-to-one relation between a set of v moments rj(si), i = 1, ..., u and the well-defined steady state response of the signal y(t) we mean 
that the moments r) uniquely determine the steady-state response of y(t). 
Let (A, B, C) be a linear system described by the equations x = Ax + Bu, y = Cx, with x(t) € R" and cr(A) S C — . Then x(t) = x p (t) + xt(f), 

with x p (t) = lim / e A ( t ~ T > Bu(t)cIt . xt(t) is the transient component of the state x(t), i.e., lim xt(t) = and x p (t) is the steady-state. 

Consequently, y v (t) = Cx p (t) is the steady-state response of the linear system (A, B, C). 



The next result shows that the models Eg and £# achieve moment matching at cr(S) and <j(Q) for all parameters G and 

H, respectively. 

Proposition 1. [^ The following statements hold. 

1. Consider the family of systems 0). Consider a v-th order model of system {!]) at <r(S) and let K(s) be its transfer 
function. Then there exists a unique G such that K{s) = B*QH(sI — S + GL)~ 1 G . 

2. Consider system (QP and the family of v-th order models (pi). Let T be the (unique) solution of equation |[3J). Then, 
for all H , any model in the family (T71) is a reduced order model of system (QP achieving moment matching at o~{Q). 

□ 

Note that the results so far cannot be applied to the case of infinite interpolation points, i.e., do not apply to s = oo. 

2.2 The MIMO case 

Consider a MIMO system of the form (fTJ), with the input u{t) £ R m , i.e., B £ C" xm and the transfer function K(s) £ 
C mxm . Let S £ C vx " and L = [h l 2 ... l v ] £ C mxi/ , k £ C m , i = l,...,v, be such that the pair (L,S) is observable. 
Let II <E C nx " be the unique solution of the Sylvester equation @. Simple computations yield that the moments 
fj(si) — K(si)li, rj(si) £ C p , i = 1, ..., v of system (fT|) at \s\, ..., Sy\ — a [Si) are in one-to-one relation with CIL Consider 
the following system 

i-Ft + a*, 

with £,(t) £ W, i]j{t) £W,G £ C" xm and H £ C pxv . The model reduction problem for MIMO systems boils down to 
finding a z^-th order model described by the equations (JSJ) which satisfies the conditions 

K(s i )l i = K{s i )l i , t=l, ...,!/, (9) 

where K(s) — H(sl — F)~ 1 G is the transfer function of §8§. The relations © are called the right tangential interpolation 
conditions, see [15] . It immediately follows that the solution to this problem is provided by a direct application of Theorem 
[21 i.e., a class of reduced order MIMO models that achieve moment matching in the sense of satisfying the tangential 
interpolation conditions Q is given by Sq = (S — GL,G,B*QTi) as in ©. 

Similarly, we may define the left tangential interpolation problem and its solution. To this end, let Q £ C ,yx ' / and 
K = [r* ... r*]* £ C xp , n £ C lxp , i = 1, ...,v, be such that the pair (Q,H) is controllable. Let T £ <C vxn be the unique 
solution of ([3]). Hence the moments r)(si) = riK(si), r)(si) £ C lxm , i = 1, ...,^, of system ([T|) at {s\, ...,s„} = o-{Q) are 
in one-to-one relation with TB. The model reduction problem boils down to finding a i/-th order model described by the 
equations |(8]) which satisfies the conditions 

riK(si) = nK(si), i — 1, ...,v. (10) 

The relations (jTUJ) are called the left tangential interpolation conditions, see [T5] and the solution to this problem is 
provided by a direct application of Theorem^ i.e., a class of reduced order MIMO models that achieve moment matching 
in the sense of satisfying the tangential interpolation conditions (TTU|) is given by E# = (Q — 1ZH, TB, H) as in ([7]). 

Note finally that also in the MIMO case the models are parameterized in L and H, respectively. Their choice is important 
in establishing appropriate directions for interpolation. Throughout the rest of the paper we discuss the SISO case, i.e., 
m = p = 1, the results being easily extended to tangential interpolation for MIMO systems. However, when necessary, we 
make specific remarks about the latter case. 

2.3 Krylov projections 

In this section we give a brief overview on the port Hamiltonian reduced order models obtained using Krylov projections. 
Proposition 2. J34. 4^ Consider the system U\). Let sq £ C be such that so $ cr((J — R)Q) and let V £ M nXL ' be such 
that V £ span{((J — R)Q — sqL)~ 1 B, ..., ((J — R)Q — sqI)~ v B}. Consider the port Hamiltonian system 

E {£=(Jr-Rr)Qr£ + BrU, 

v '\ib = B;.Q r t 

with 

J r = V*QJQV, R r = V*QRQV, 

(12) 
Q r = (V*QV)- 1 , B r = V*QB. 

Then £y matches the first v moments of $7$) at sq. □ 



Furthermore, in [23] a Krylov projection technique is proposed for the computation of a port Hamiltonian reduced order 
model that achieves moment matching, i.e., tangential interpolation, in the MIMO case. 

Theorem 3. [23] Consider a MIMO system Q) with the input u(t) £ W m and the output y(t) £ K m , i.e. B £ C nxm . Let 
si, ..., s u £ C be such that s\ =/= 8% =/= ... =/= s u . Let l\, ..., l v £ C™ 1 and 

V = [{ Sl I - (J - Q)R)- l Bh ... (s u I - (J - Q)R)- l BQ. (13) 

Let M £ C XI ' be any nonsingular matrix such that V — VM is real. Define W — QV(V*QV)~ 1 and 

J r = W*JW, R r = W*RW, 

~ J ^ (14) 

Qr = (V*QV)- 1 , B r = W*B. 

Then the reduced order model So as in ill]) is port Hamiltonian, passive and satisfies the right tangential interpolation 
conditions (0). □ 

Note that Theorem [3] yields a class of reduced order port Hamiltonian models, with state-space realizations parameterized 
in M, which achieve moment matching. Results which extend this technique to matching a number of moments which is 
twice the number of selected interpolation points can be found in [41] . 

Remark 1. To the best of our knowledge, the existing moment matching based model order reduction techniques do 
not yield approximation error bounds. However, a choice of interpolation points that improves the accuracy of the 
approximation, i.e., stems from solving an H2 optimal approximation problem, is at the mirror images of the poles of the 
approximant. In detail, if K(s) is the approximation that yields the best approximation of K(s) associated to ([T|), in the 
B.-2 norm, then K(s) satisfies the conditions 

K(-Xi) = K(-Xi), (15a) 



dK{s) 



ds 



dK(s) 



ds 



(15b) 

s=-\i 



with i = 1, ..., i/, where Xi are the poles of K(s), see e.g., [331 I3S1 I2D]- Hence, an accurate approximation of a system ([T]), 
in the sense that the T^-norm of the approximation error decreases, is obtained by finding a reduced order model that 
achieves moment matching in the sense of relations (|15[) . Since the poles of the approximation are not known in advance, 
a suitable choice of interpolation points (see [20J [21]) is 

Si = -Aj, i = l,...,v, (16) 

where Xi £ a ((J — R)Q), with highest residues of the transfer function K(s) . D 

Proposition [JJ directly applies to Markov parameter matching, i.e., matching at infinity. However, this is not the case 
for Theorem [3j This is due to the fact that the construction and direct application of the projectors W and V does not 
yield a reduced order model with port Hamiltonian structure. However, in [401 142] a solution is obtained in a specific 
set of coordinates, i.e., performing a coordinate transformation the port Hamiltonian system is brought in a form where 
Proposition [3] applies. A precise definition of Markov parameters and the corresponding detailed arguments are found in 
Section [3J 

To the best of the authors' knowledge there is no structured way of determining the number v of interpolation points 
required for approximation. However, a fair choice appears to be the number of large Hankel singular values of the given 
system, see, e.g., [3"0"] . 

3 Markov parameter matching as moment matching for a class of descrip- 
tor systems 

As mentioned in Section 12.11 the time-domain moment matching tools do not directly apply to s — 00. In this section 
we present a time-domain approach to moment matching at s = 00, i.e., Markov parameter matching. The idea is to 
perform the change of variable r = 1/s O s = 1/r and turn the problem of Markov parameter matching into a moment 
matching problem at r = 0. First we discuss the problem of moment matching at any r £ C and provide time-domain 
interpretations of the notion of moment and moment matching. Based on this notion of moment, we characterize the 
classes of parameterized reduced order models that match the moments of a given system at r and in particular at r = 0. 

Consider a linear, minimal system described by the equations 

X = Ax + Bu, 

r ( 17 ) 

y = ex, 



with x(t) £ K", y(t) G M, u(t) G M. Let K(s) = C(s/- A)" 1 ^ be its transfer function. The first v+1 Markov parameters 
arc the coefficients of the series expansion of K(s) around s — oo, i.e., they are the first v + 1 moments of K(s) at oo, 
namely 

r?o(oo) = 0, r] k (oo) = CA k ~ 1 B, k = 1, ..., v. (18) 

Let r € C and note that the matrix pencil I — At is regular, i.e., det(7 — At) + 1 0, for some r (see e.g., [ID])- According 
to, e.g., [5J> (I — ^4r) _1 exists and so the function K(t) = K (~) = C(I — At)~ x Bt is well-defined. Furthermore, we have 



that 



jfc+i 



dr /c+i 



" ' ' ' (fc + 1)!C[(J - ^r)- fe - 1 J 4 fc B + (7 - Ar)- fc - 2 A fc + 1 Br], yielding 

P = C(/-Ar)- fe - 2 A fe S. 



1 d fc+1 7s:(r) 



(fc + 1)! rfr fc + 
The moments of K(t) at t = t* G C are given by 

%(r*) = 



1 



d k+1 K(i 



and the moments ?7o(oo), ...,?7„(oo) are given by 



77fc(oo) = 



(fc + 1)! dr fe +! 



1 d fe+1 if(r) 



(19) 



(fc + 1)! dT k + x 



= rjk(0). 



T = 



We now consider the following moment matching problem. Given the function K(r) and the point t* G C find K(t) such 



that the first v + 1 moments at r* match, i.e. ^(r) = dT k+i ( r *) = ^fit(t*), for all fc = 0, ..., ia In particular, we are 
interested in the case r* = 0, which solves the Markov parameter matching problem. 



Proposition 3. Consider the system |J7| ) and r* G C Let 

L=[l ... 0] eR lx (" +1 ', 



5' 



T* 1 

t* 1 



r* 1 
. . . t* 



c (H-i)x (i/+i)_ 



Assume that At* 7^ 1, /or any A G cr(A). 

i. Let LI G C" x( '''' +1 - ) &e the unique solution of the Sylvester equation 

AUS + BL = U. 
Then the moments T)o(t*), . . . , rj v {T*) of K(t) = K(l/r) satisfy 

[?} (r*) Vi(r*) ... rju(r*)} = CTLS. 

2. Let LT G C nx ^ +1 ' be the unique solution of the Sylvester equation 

AUS + BLS = n. 
Then the moments 5Jo(t*), . . . , r]v{T*) of K(r) — K(1/t) satisfy 

[rj (r*) rjiir*) ■•■ ^(r*)] = Cfi. 



(20) 



(21) 



□ 



Proof. The assumption that At* 7^ 1, for any A G cr(yl) and the regularity of the matrix pencil L — At yield that the 
solutions LI and LT of the Sylvester equations ([2"0"]) and (j2"Tj) . respectively, exist and are unique (see e.g., [TO]). 
To prove the statement (1), let U — [U a LIi . . . n„]. By (f2"0"|) we have the sequence of equalities 

An T* + b = n <=> n = (I - At*)- 1 ^, 
ail, + AniT* = n x «*■ n x = (/ - At*)- 1 ^ 

= (7 - At*)- 1 ^^ - At*)~ 1 B 
= {L-At*)- 2 AB 1 



An,_! + ail,t* = n„ <*■ n„ = (I - At*)- 1 !^ 

= (/-AT*)-' / - 1 yl l/ B, 



from which the claim follows. To prove the second statement, let II = [IIo II i ... TL V ]. By (f^T|) we have the sequence of 
equalities 

ATL t* + Bt* = S «*• fl = (I- Ar*) _1 Br*, 
Ati + Ati lT * + B = Hi ^H, = (J-Ar*) _1 (Ano + B) 



An^_i + Afl^r* = n„ *>- n v = (j - At*)- (u+1) a u Bt* 

+ {I-At*)- u A u - 1 B, 



from which the claim follows. 
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Proposition [3] can be extended to the general case of any (non-derogatorjl 2 ) matrix S that satisfy the assumption made on 
cr(A). Hence, throughout the rest of the section we make the following standing assumption. 
Assumption 1. A/i 7^ 1, for any A £ cr(A) and /1 £ o~(S). 
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Figure 1: Block diagrams of systems (J22D (a) and ([23]) (b). 



Example 1. Consider the system depicted in Fig. 1(a) described by the equations 



xi = x 2 , 

±2 = —axi — bx2 + u, 

y = xi, 



(22) 



with the transfer function K(s) = ' 



s 2 +bs+a ' 



Consider now the system depicted in Fig. 1(b) where we have replaced the integrators with differentiators, i.e., t = -. 
Note that x\(t) = 4f(%2{t)) and X2(t) = 4z(—axi(t) — bx2{t) + u(t)) yield the descriptor state-space realization 



x 2 = xi, 
-aii — 0x2 = Xi — u, 
V = *ij 



(23) 



with the transfer function K(t) = K (^j = a7 . 2 \ T+l - Hence, the Markov parameters of (|22|) are the moments of the 
system (QSD at r = 0. D 







- u 


Ax — x - 
y = Cx 


-Bit 




u — Slj 
8 = Lu 


6 = 


y 















Figure 2: Interconnection of the signal generator and the system (|2l)]) . with transfer function K{t). 

The following result provides a time-domain interpretation of the moments of K(t) = K(1/t), where K is the transfer 
function of the system (fTT)) . Namely, the moments of K(r) are related to the well-defined steady-state response of the 
output of a descriptor state-space representation associated to K(t), with the input given by a signal generator defined 
by the interpolation points given by the set a(S) (see Fig. [2J. 
Theorem 4. Consider the system J_?7|) and let S be such that o~(S) C C . Assume that o~(A) c C~. Let 



lo = Suj, w(0) ^ 0, 



(24) 



withoj{t) G W. 



A matrix is non- derogatory if its characteristic and minimal polynomials coincide. 



Consider the descriptor state-space representation associated to the transfer function K given by 

Ax = x — Bu, 
V = Cx, 



(25) 



with x(t) £ l n , u(t) £ R and y(t) £ R. Then the moments rjk(Ti), k = 0, 1, 2, ..., i — 0, ..., v , of the system \llty at 
o~(S), are in one-to-one relation with the well-defined steady-state response of the output of the interconnection of 
systems \25}) and {21$, with u — Lu>, where L is such that the pair (L, S) is observable. 



(26) 



ii. Consider the descriptor state-space representation associated to the transfer function K given by 

Ax = x — Bu, 

y = Cx, 



with x(t) £ R n , u(t) £ R and y(t) £ R. Then the moments rjk(Ti), k = 0, 1, 2, ..., i = 0, ..., v , of the system |_/7[ ) at 
c(S), are in one-to-one relation with the well-defined steady-state response of the output of the interconnection of 
systems \2b}) and \2J$, with u = Luj, where L is such that the pair (L, S) is observable. □ 



Proof. By the interconnection between ([2"5)l and (pM|) . the first equation in ([23)) yields x — Ax + BLu = Ax + (II — AHS)lo, 
where II is the unique solution of ([2"0")) . This further yields x — ITw = A(x — IIcj). The output y(t) satisfies y(t) = 
C(x — ITcj) + CUdi. The assumptions <j{A) c C _ and cr{S) C C° yield lim [x(t) — IIu;(£)] = 0, implying that the first term 

t— ^oo 

of y(t) vanishes. Hence the steady-state response exists and is well-defined and equal to Clio; = CUSui, which proves the 
claim i. Since cr(A) c C~, the statement ii. is a direct application of Theoremjl] with the signal u(t) replaced by u(t). D 

Corollary 1. Consider the system Ji7| ) satisfying the assumptions in Theorem^A Then the first v Markov parameters of 
|17[ ) are in a one-to-one relationship with the well-defined steady-state response of the system i25\) (or h2b}) ) to the input 
given by a signal generator defined by the interpolation points t* = with multiplicity v . □ 

We are now ready to define the notion of reduced order models that match the moments of K{f) — K(1/t). 
Consider an observable pair (L, S) and the system 



-n 



f = F£ + Gu, 



(27) 



with £(i) £ W and II = II, or II = II. Let K n {s) be the transfer function of system ([27]) and let K n (r) = K n {l/r). 
Then, if II = LT and A/i ^ 1, for any A £ a{F) and /i £ o-(S), (|27l) is a reduced order model that matches the first v 
moments of K(t) at cr(S) if there exists an invertible matrix P £ R" xjy such that 



CnS 1 = HPS, FPS + GL = P. 
If t* = 0, then the Markov parameters of Ku(t) match the first v Markov parameters of K(s). 



(28) 



Finally, for II = IT, the system ([2"T[) is a reduced order model that matches the first v moments of K(t) at cr(S) if there 
exists an invertible matrix P £ M VX1 ' such that 



Cn = HP, FPS + GLS = P. 



(29) 



Remark 2. Assume S is invertible, i.e. r* ^ 0. Let P — I. A reduced order model that matches the moments of K (r) 

at t* is given by equations ([2T|) with F = (I — GVjS^ 1 and H = CIL Furthermore, if P = I, according to ([2U)) another 

reduced order model that matches the moments of K(t) at r* is given by equations ([2"7]) with F = S^ 1 — GL and H = CIL. 

D 

Example 2. Let 

"010 

L= [1 0], S= 1 





CII= [771 r? 2 r? 3 ]. 



The first three Markov parameters of K(s) are CITS' = [0 771 772] = [0 CB CAB]. A reduced order model that matches 
these Markov parameters is given by equations ([2"5t , with 



F = 





M 
/31 



G = [10 0]*, H = [771 % 773]. 

In general, let lxv = [0 ... 0] £ R lxu , Si = [1 ... 0] £ R lxv and 



1 



1 /33 



s 2 = 



T ' 





T 





Then L = [1 Oi XI /] and S 



Si 
0, S 2 



are such that the pair (L, S) is observable. Furthermore, let CII — \r\\r\2 ■ ■ ■ f]v Wv+i] € 



l lx '" +1 '. The first ^ + 1 Markov parameters are 

GILS' = [0 J?! .. .nA = [° c^ c- 45 ■•• c^^s]. 

Let F = n , G = [Gi G2]* and P = I. Solving (f^Hj) yields a family of reduced order models of dimension 

-t*21 -f*22 J 

;/ + 1 that achieve Markov parameters matching characterized by G\ = 1, G2 = 0, Fn = 0, F12 — [0 ... 1], F12, 

iJ = [vi V2 ■ ■ ■ Vl> Vu+i] an d F22 free parameters. 

Finally, let ft = [ft ft], where ft = [fti ... Iljel"". In this case, © becomes 



[0, AU] 



Si 
0, 5 2 



[5 n> 



Si 
0, 5 2 



= Po n], 



yielding the equations 



ft = 0, AUS 2 + BSx = n. 



(30) 



Note that the first Markov parameter is 770 (00) = Gflo = and the remaining v Markov parameters are in a one-to-one 
relation with Gil. A reduced order model of dimension v that matches the first v + 1 Markov parameters, i.e. 770(00) = 
and 771(00), . . . , ?7„(oo), is given by system (j2"T)l if there exists an invertible matrix P G E jyx,y such that 



CU = HP, FPS 2 + GSi = P. 



(31) 

□ 



3.1 Alternative approach 

In this subsection we give an alternative interpretation of the notion of moment matching for K, in the sense of the 
approach taken in [7]. The moments of K(t) are in one-to-one relation with the well-defined steady-state response of the 
asymptotically stable descriptor systems (f^5|) or (|2T)1) driving a generalized signal generator. This approach yields new 
families of reduced order models that achieve moment matching at r and in particular, Markov parameter matching when 
r = 0. 
Proposition 4. Consider the system JJ7[ ) and r* £ C. Let 



Q 



T* 








1 


T* 










1 















T* 



(j/+l)x(v+l) 



be such that Assumption^ holds and 



ft= [1 ... 0]* £ 



lU+l 



1. Let T G (^(f+ijx™ fr e 7/jg unique solution of the Sylvester equation 

T = QTA + KC. 
Then the moments of K(t) = K(1/t) at r* satisfy 

[Vo(r*)m(r*) ... rj v (r*)]* = QTB. 

2. Let T G (^(f+ijx™ fr e 7/jg U nique solution of the Sylvester equation 

T = QTA + TICA. 
Then the moments of K(t) = K(1/t) at r* satisfy 

Mr*) Vi(r*) ... rj»(r*)}* = Q(QT + TZC)B. 



(32) 



(33) 



□ 
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Proof. The assumption that At* ^ 1, for any A <E cr(A) and the regularity of the matrix pencil / — At yield that the 
solutions II and fl of the generalized Sylvester equations ([521) and ([3"3")l . respectively, exist and are unique (see e.g., [TO]). 
To prove the statement (1) let T = [Yj$ T* . . . Y*]*. From (l32j) . we have the sequence of equalities 

r*T A + C = T ^T = C(I- At*)' 1 , 
T Q A + t*T x A = Ti «*• Ti = T Q A{I - At*)- 1 

= C(I - Ar*)- 1 ^! - At*)- 1 
= CA(I-At*)- 2 7 



n„_iA + t*T,A = T v *> T v = T w _i(J - At*)- 1 

= CA v (I-At*)- u - 1 , 

from which the claim follows. 

To prove the second statement let T = [T*, T* . . . Y*]*. From ([55)1 we have the sequence of equalities 

r*¥ Q A + CA = %&T = CA{I - At*)- 1 , 
f Q A + t*T x A = % <& Ti = Y A(I - At*)- 1 

= t*CA 2 (I-At*)- 2 , 



T„_iA + t*T u A = T V ^T V = t*CA v (I - At*)-", 



from which the claim follows. 



□ 



Remark 3. Let Q and TZ be any two matrices such that the pair (Q,7Z) is controllable. Then the moments described by 
([2]) are parameterized by the elements of TZ. However, there exists a coordinate transformation T such that T- 1 QT = Q, 
as in (gj) and T^K = TZ, as in (gj, yielding the moments of ^T(t) at ct(Q) = a{Q). D 



j4i = x — Bw 
j/ = Cx 



v = y 



Co = Quo + IZv 



Tx 



Figure 3: Interconnection of the signal generator and the system (1261) . with the transfer function K(t). 

Consider the generalized signal generator 

Co = Quo + TZv, w(0) =0, (34) 

with uo(t) € W, v(t) G K. Consider the interconnection v = y (as in Fig. [3]) between the signal generator and the 



system (f2"6")l . Let d = cj — Tx. The signal d(t) satisfies d = Co — Tx =£■ d = Qw + (TvLC — T)i. This further yields 
d = Quo — QTAx = Quo — QT(x — Bu). In conclusion d satisfies the equation 

d = Qd + QTBu. (35) 

Conversely, let T be such that d — uj — Tx satisfies equation (|35p . Then T satisfies Qd+QTBu = Co — Tx — Quo+lZCx—Tx, 
which yields QTx — QTBu = (JZC — T)x. Hence QTAx — (1ZC — T)x, for all x, i.e., T satisfies equation (j32|) . 
Consider the generalized signal generator ([M]) interconnected through v — y (as in Fig. 2]) with the system (|2"5")) . Let 



Ax = x — Bii 
y = Cx 



v = y 



Co = Quo + IZv 



d = Q(oj - Tx) 



Figure 4: Interconnection of the signal generator and the system (|2"5"1) . with the transfer function K(t). 

d= Q{uo-Tx). The signal d{t) satisfies d= Q{Co-Tx) = Q(Quo+ (RCA- T )x + KCBu) = Q(Quo- QT Ax + KCBii) = 
S(Sw - QT(x - Bu) + IZCBii), hence 

d=Qd+Q(QT + KC)Bu. (36) 

Theorem 5. Consider system J_?7| ) and Zet Q 6e smc/i i/iaf cr(Q) C C . Assume that (J (A) C C~ and x(0) — 0. Let the 
system J5*^[ ) 6e smc/i f/iai £/ie pair (Q,7Z) is controllable. 

i. Consider the interconnection between the system V25\) and the system \3$ , defined by v — y, with the output d(t). 
Then the moments rjk(ji), k = 0, 1, 2, ..., i = 0, ..., f, o/ i/ie system {lty at o{Q), are in one-to-one relation with the 
well-defined steady-state response of the output d(t), for u = S(t). 
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ii. Consider the interconnection between the system A26}) and the system \3J$ , defined by v — y, with the output d(t). 
Then the moments rjk(ji), k = 0, 1, 2, ..., i = 0, ..., v, of the system J-?7| ) at cr(Q), are in one-to-one relation with the 
well-defined steady-state response of the output d{f), for u = \(t), where \{t) denotes the Heaviside step function. □ 

Proof. The proof is similar to the proof of Theorem [T] (see also [7]), hence, we provide a sketch for the statement i. By 
the assumptions, the steady-state response of d(t) is well-defined. By equation (|35j) . since u(t) = 5(t), we have that 
D{t) = (tI — Q)~ 1 QYB, where D(t) denotes the Laplace transform of d(t). D(t) contains the moments in the terms 
K(t*)H/(t-t*), where t* € a(Q) and K(t) is the transfer function of the system $ZJ$. Recalling that i(t) = $$■ = 5(t), 
the proof of the statement ii. is identical to the proof of statement i. □ 

Corollary 2. Consider the system fiTf ) satisfying the assumptions in Theorem\5\ Let o~(Q) = {0, ..., 0}. Then the first v 
Markov parameters of \17\j are in a one-to-one relation with the well-defined steady-state responses of the signals d(t) and 
d(t), to inputs u(t) = 5(t) and u(t) = l(t), respectively. □ 

We are now ready to present new families of reduced order models that achieve moment matching at c(Q). 
Consider a system described by equations of the form 

with f(f) e W, u{t) £ R, ijj(t) e R, and E, F, G, H of appropriate dimensions. Let K T = H(sE - F) _1 G be the 
well-defined transfer function of (|3"T)) and let K*c(t) — K-c{1/t). 

Let the signal C,(t) be such that ( = uj — PE^. Then, the moments of the transfer function K~r(r) of (|57|) match the 
moments of the transfer function K(t) of flTT|) at a(Q), in the sense of Theorem^ if £(£) satisfies C = QC, + QTBu, which 
is equivalent to the existence of an invertible matrix P such that 

QPF + KH = PE, QPG= QTB. (38) 

Similarly, let £(£) be such that C, = Q(w — PE£). Then, the moments of the transfer function K~c(t) of (|37l) match the 

moments of the transfer function K(t) of (fTT)) at er(Q), in the sense of Theorem[SJ if £(£) satisfies £ = QC+2(2"f +7^C)-E? U : 
which is equivalent to the existence of an invertible matrix P such that 

QPF + KH = PE, Q(QP + KH)G = Q(QT + KC)B. (39) 

In the sequel, we present particular instances of E, F, G and H, which satisfy the matching relations (|38l) and (|39"|) . 

respectively. 

Proposition 5. Consider system \11ty with the transfer function K(s). Let K(t) = K(1/t) and let T and T be the unique 

solutions of the Sylvester equations A32\) and (0), respectively. Furthermore, let Q and 1Z be such that the pair (Q,1Z) is 

controllable. Then the following statements hold. 

1. A reduced order model with transfer function Kx(t) which matches the moments of K{t) at o~(Q) is given by the 
system (37|) with E = Q - KH, F = L and G = -TB. 



2. A reduced order model with transfer function Kx(t) which matches the moments of K{t) at o~(Q) is given by the 
system (37|) with E = Q- KH, F = L and G = -QTB and output ip = H£. 

3. A reduced order model with transfer function Ky(r) which matches the moments of K(t) at o~(Q) is given by the 
system fSTjf with E = Q- KH, F = L and G = -(QT + KC)B. 

4- A reduced order model with transfer function K~r(r) which matches the moments of K(t) at o~{Q) is given by the 
system (37|) with E = Q- KH, F = L and G = -Q(QT + KC)B and input ii(t). 

Proof. We only prove statement (1), since the statements (2)-(4) follow similar arguments. Let K(t) = C(I — At)~ 1 Bt 
andletK Y (r) = H(tI-Q+KH)- 1 TBt. Hence E(t) = K(t)-K t {t) = [C-iJ(T/-Q+^iy)- 1 T(/-^r)](/-AT)- 1 Br. 
Exploiting equation (fTJ) and performing some algebraic computations yield E(t) — E(t)(C + HTA)(L — At)~ 1 Bt, where 
E : C -> C, E(t) = 1 - H{tL - Q + KH)- x K. Note that 

%) = 1 H(rL-Q^K 



l + HirL-Q)-^ l + Hirl-Q^K' 

It follows that for all r,; £ cr(Q), i = 1, ..., v, Eirf) — and then E{ri) — 0, which proves the claim. The rest of the claims 
are proven in a similar way. □ 

Remark 4. Let F = I. Then a reduced order model that achieves moment matching of K{t) at <r(Q) as in Proposition 
is obtained for P = —I or P = —I, where P uniquely satisfies (|55)) and P uniquely satisfies (|39p . respectively. Note 
that the uniqueness of the solutions holds if and only if a(Q) n u(Q — KH) =0. □ 
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Example 3. Let 



Q 





1 
1 



n 



and consider a reduced order model (|37[) with E = I , G = [g± gi 53]* and H = [hi h^ hs]. Let T be the unique solution of 
J]} and let TB — [770 r/i 7/2]*. System (l37l) matches the first three Markov parameters of the system (fl7|) . i.e., the moments 
QTB = [0 770 171]*, for any P that uniquely satisfies (|3"5j) . Choosing P = I, and solving (f3"5f yields 



with the transfer function 



F = 






1 







% 








1 


, G = 


Vi 


-a 


-6 


— c 




93 



H = [l 0], 



«"(*) 



i] s 2 + (771 - ?7oc)s - 7706 - ryic + .93 



s 3 _|_ g 2 c _|_ S J _|_ fl 

where a, 6, c and 53 G R are free parameters. In general, let 



Q = 



I, 



(u-l)x(u-l) 



i) 



, K = 







and E = I. Choosing P = I and solving (1381) yield a class of reduced order models that match v Markov parameters of 
the system (fT7|) described by 

" /„_! 



F 



/1 



F,* 



, G 



.9 



, if = [1 



% 



with F 2 G W- 1 and /1 G R and g G R free parameters, and S such that QTS = [0 3*]*. Note that A/J = 7^ 1 for all 
A G c(F), since = n £ ct(<2), and thus P = I uniquely satisfies (|38p . D 

Example 4. Let Q G R l/XI/ and H G R 1 ' be such that the pair (Q,1Z) is controllable and let T be the unique solution of 
equation ([TJ) . Furthermore, let P — I be the unique solution of ([55)1 and assume E = I. Then a class of reduced order 
models that match the moments of K(t) at a(Q) is given by 

F = (Q + KH)-\ G = (Q + UHy^QT + KC)B, 

with H a free parameter, such that Q + 1ZH is invertible and A/i 7^ 1, for any A G a(F) and /j G c(Q). 

In particular, let Q be a Jordan block of order 1/ with all eigenvalues at zero, 7Z — [1 ... 0]* G R", (QT + 1ZC)B = 
[rjo 771 ... 7?i/-i]* and if = [/ii ft-2 ••• ^1/] G R lx,y . Then det(Q + 1ZH) = h v . Choosing h v 7^ 0, we obtain a class of reduced 
order models described by equations (|37|) and parameterized by h\, ..., h v , that matches the first v Markov parameters of 
the system (fTTJ), with 



F 






1 





" 












1 





, G = 


Vi 











1 




Vu-1 


1 


ha 


h-2 


h v J 




I 9 



and g = 



Va 



m=i h iVi 



h v 



Note that by construction $■ ct(F) and so A/i = for all A G cr(F), since = /1 G c(Q), i.e., 



P = I uniquely satisfies (|3U)) . 



D 



3.2 Relation with Krylov methods 

In this section we establish a connection between the Krylov projection in the Arnoldi method (see, e.g., [2j|42]) and the 
solution of the Sylvester equation (|20p . We will only sketch the result which follows arguments from the proof of Lemma 
lETTl Consider the system ([17} and let V G C" XI/ and W G C nKU be such that 



V 

W 



[B AB ... A V - 1 B], 

[C* {CA)* ... (CA"- 1 )*} 



(40a) 
(40b) 



Let S be as in Proposition [3] and note that AVS + BL = V, with L such that the pair (L, S) is observable. Let T be an 
invertible matrix such that S = TST~ l , where S G C UXI/ is any matrix such that a(S) = a(S). Hence AVTS+BLT = VT 
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yielding that the unique solution II of (OH)) satisfies II = VT. Similarly, II = VT, where II is the unique solution of (|2ip. 
Note that similar results are obtained for T, the unique solution of (J22) or ([55]) . i.e., T = TW*, where T is an invertiblc 
matrix such that Q = TQT^ 1 , with Q G C yxu as in Proposition ffl 

The following result establishes a relation between the class of z^-th order models that match a set of prescribed Markov 
parameters, obtained using Krylov projections and the class of reduced models En- A similar result is obtained for a class 
of models E-r- 
Proposition 6. Consider the system {lty and the system 

£ = W*AV£ + W*Bu, 

with £(£) G W . The following statements hold. 

1. Let V be as in \40aty and let W be such that W*V—I. Let Ejy be the family of v order models of J-/7] ), described by 
equations 071) and parameterized in W . Then, £vy = Sri, with Sn as in \21ty , with II G {II, n}, where II is the 
unique solution of i20\) and n is the unique solution of \21}). Analogously, let W be as in &40b\ ) and let V be such 
that W*V—I. Let Ey be the family of v order models of |_?9f ), described by equations 071) and parameterized in V. 
Then, Ey = Sn- 

2. Let W be as in RJJM) and let V be such that W*V=I. Let Ey be the family of v order models of |J7[ ), described by 
equations 071) and parameterized in V . Then, Ey = Ex, with Ex as in \31ty , with T the unique solution of {32]) . 
Analogously, let V be as in J40a\ ) and let W be such that W*V—I. Let E^y be the family of v order models of J_?7| ), 
described by equations OT1) and parameterized in W . Then, Y<\y = Ex- D 



Proof. Proof of Statement (1). Consider a model Y,w and let II = VP, with II the unique solution of (P2U1) and P an 
invertible matrix. Thus, CVP = GIL which yields CIIS = CVPS. Furthermore, note that W*AVPS + W*BL = 
W*{AILS + BL) = W*VP = P for any W such that W*V = I. Hence, the conditions ([55]) are satisfied and the claim is 
proven. Similar arguments hold for II = n, the unique solution of (|2lT) . as well as for the case Ey. 



Proof of Statement (2). Let T be the unique solution of (|55|) and consider the family of v-th order models Ex, as in (|57|) . 
Without loss of generality assume E — L in (|57|) . Consider a model Ey and let T = PW* with P an invertible matrix. 
Hence, QPW*B = QTB. Furthermore, note that QPW* AV + TZCV = (QTA + TZC)V = TV = PW*V = P, for any V 
such that W*V = L. Thus, the conditions (1551) are satisfied and the claim is proven. Similar arguments hold for the case 
En/. □ 

Note that the results from Proposition [5] do not apply to the family of models E^, with T the unique solution of (|55j) . 

4 Matching with preservation of the port Hamiltonian structure 

4.1 Matching at finite interpolation points 

In this section, we solve the following particular instance of theQ] at a prescribed set of finite interpolation points. 
Problem 2. Given a linear, port Hamiltonian, SISO system (fT]), find the observable pair (L,S), where L s C lx " and 

5 G ( C V/V , such that a(S) D o~((J ~ R)Q) = and the free parameter G G C 1 ', such that the following properties hold. 

(pi) The family of systems (S — GL, G, B*QIT) parameterized in S, L, G, described by equations ([5]), with n the unique 
solution of @, match the moments of (JTJ) at a(S). 

(p2) There exists G such that the family of systems (S — GL, G, B*QH) parameterized in S and L, described by equations 
(|r?j) preserve the port Hamiltonian structure. 

(p3) The family of systems (S — GL,G,B*QH) parameterized in S, G, L, described by equations ([6]) are accurate 
approximations of (JT|), in the i^-norm sense. 

(p4) The computation of the family of systems (S — GL,G,B*QH) parameterized in S, G, L, described by equations 
([5]) is calculated efficiently, i.e., avoiding the computation of the moments and the solution n of @. □ 

The solution to subproblcm (pi) is provided by Theorem [2] which characterizes the family of systems (S — GL, G, B*QH) 
parameterized in S, L, G. Our main goal is to solve subproblem (p2). In short, fixing S and L we compute the is-th 
order system (S — GL, G, B*QU) that preserves the port Hamiltonian structure and matches the moments of ([T|) at a(S). 
This model characterizes a class of state-space models parameterized in S and L, which have the same transfer function. 
The parameters S and L can be selected to satisfy further constraints such as additional physical structure and improved 
accuracy of the approximations, by solving subproblem (p2) at the mirror images of a set of poles of the given system (fT]), 
i.e., the port Hamiltonian reduced order model satisfies conditions (I15a|) . thus shortly addressing subproblem (p3). In the 
case of MIMO systems, a proper selection of L defines the desired directions for the solution of the tangential interpolation 
problem, i.e., the reduced order port Hamiltonian system satisfies the right /left tangential interpolation conditions © or 
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pop , respectively. Although it is not the primary goal of this paper, subproblem (p4) is addressed using the result in 
Theorem [3] and the relation between the projections and the solutions of the Sylvester equations @ and Q, provided 
by Lemma TB. II Hence, the output B*QU of (|6]) is computed using efficient numerical algorithms stemming from Krylov 
subspace techniques as in, e.g., [25] , 

Note that a similar problem and results are obtained using the family of systems (Q — 1ZH,TB,H) parameterized in 
Q,1Z, H, as in 0, with T the unique solution of Q. 

Our main focus is solving subproblem (p2) of Problem[2] The results consist of families of reduced order, port Hamiltonian 
models which are subclasses of the class of models Eg in @ and E# in (J7J|, respectively. The port Hamiltonian models 
match the moments of the given port Hamiltonian system and possess parameterized state-space realizations. Note that 
all the model from the subclass have the same transfer function. The parameters can be used for enforcing additional 
structure such as, e.g., specific Hamiltonian functions and/or diagonalized dissipation. In the MIMO case, the parameters 
can be used to find the reduced order models that satisfy the tangential interpolation conditions. 

Consider the linear, SISO, port Hamiltonian system (TTJ). Let (L, S) be an observable pair and let (Q,1Z) be a controllable 
pair. Consider the families of reduced order models Eg as in © and E# as in ([7]), respectively. Throughout the rest of 
this section we make the following working assumption. 

Assumption 2. The matrix Q is invertible. Furthermore, a(S)r\a((J-R)Q) = 0, cr(Q)na(A) = 0, a(S)r\a(S-GL) = 
and c(<2) n a(Q — 1ZH) = 0, i.e., the interpolation points are not among the poles of either the given system or its 
approximations . 

Note that by Assumption^ the families of models Eg and E# are well-defined. 
Proposition 7. Consider the port Hamiltonian reduced order model given by 



£ = (J - R)Q£ + Bu, 
ijj = B*Q£, 
with £(£) £ W . Then, the following statements hold. 

1. Let n be the unique solution of equation |p|) and define 

J=n*QJQn, R = TTQRQU, 

Q = (n*gn)- 1 , b = itqb. 



(42) 



(43) 



Let En be a port Hamiltonian model described by equations H^l and G3I). If o~(S) D cr((J — R)Q) = 0, then En 
matches the moments of the system (QP at o~(S). 



2. Let T be the unique solution of equation |2P and define 

J = TJT*, R = TRT*, 

Q = (TQ -1 T*) - \ B = TB. 



(44) 

Let Ex be a port Hamiltonian model described by equations CM) and \44\ l- If °~(Q) ri cr((J — R)Q) = 0, then Ex 
matches the moments of the system (QP at cr(Q). □ 

Proof. Let K(s) be the transfer function of system ([!]). 

Proof of Statement (1). Let Kjj(s) the transfer function of the system En- Assuming that o-(S) D cr((J — R)Q) = 0, the 
result is obtained by simply checking the moments at Si € o~(S), i = 1, ..., v, i.e., checking that K(si) = Ku(si). To this 
end, note that K u (s l ) = B*Q(s t I ~ (J - RjQ)- 1 !! = B*QU(U*QU)- 1 (s l I ~ U*Q(J - R)QU(U*QU)- 1 )- 1 n*QB. Since 
n is the solution of ©, (J- R)QU = IIS-BL and hence K n (s l ) = B*QIl(sJ - S + (WQn^QBL^iWQlVj^QBL. 
By Proposition [TJ for A = (n*<5n) _1 Q_B, we have that K(si) = Kn(si), which yields the result. 

Proof of Statement (2). Let K T (s) be the transfer function of the system E T . Note that K T {s t ) = B*T*(TQ- 1 T*)- 1 (s i I- 
T(J - i?)T*(TQ- 1 T*)- 1 )- 1 TB. Since T is the solution of ©, we have T(J - R)Q = QT - TZB*Q and hence K(s t ) = 
B*T*(TQ- l T*)- 1 (s i I -Q + TZB^^TQ^T^-^-^B. By Proposition HJ for H = B^^TQ^T*)- 1 , we have that 
K(si) — Ky(si), which proves the result. □ 

Remark 5. Let (|42|) be a reduced order model of |T]). Then, by Theorem [I] the model (|42|) matches the moments of (flj 
at a(S) if 

B*Qn = B*QP, (45) 

where F £ l" x " is an invertible matrix such that 

(J-R)QP + BL = PS. (46) 

Note that equation (|4"5|) is satisfied by 

P = Q- 1 = n*QU. (47) 

Since we assume that a(S) D a(( J — R)Q) — 0, then P as in (|4"T|) is the unique solution of the Sylvester equation (|4"6")l (see 
also 0). □ 
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Remark 6. Let systems En as in (T4"2"j) and Ey as in (fTTj) be two reduced order port Hamiltonian models that match the 
moments of |T]), respectively. Then, by Lemma IB. II they are equivalent^, i.e., there exists an invertible matrix T such 
that UT = V and J r = T*JT, R r = T*RT, Q r = T-*QT~ X and B r = T*B, hence Ey = E n r and furthermore, Ey is a 
member of the class of models Eg . □ 

Note that the result in Proposition [7] shows that there exists a port Hamiltonian model that matches the moments of a 
given port Hamiltonian system by direct computation of the matrix n. However this can be avoided by showing that the 
model En is a member of the class of reduced order models Eg that match the moments of ([T]) at cr(S), for a particular 
instance of the parameter G. Similarly, we also show that the model Ex is a member of the class of reduced order models 
T,h that match the moments of (JTJ) at <r(Q), for a particular selection of H. 
Theorem 6. Let Assumption^ hold. Then the following statements hold. 

1. Let Eg, as in ^\), be a reduced order model of the port Hamiltonian system {!]], matching the moments at o~(S). 
Then Eg is equivalent to a port Hamiltonian system En, as in \^2fy and \4Sfy , i.e., Eg preserves the port Hamiltonian 
structure of the system (T7J), if and only if G — (XI* QTl)~ l Tl* QB . 

2. Let Tih, as in (p|), be a reduced order model of the port Hamiltonian system ([IP, matching the moments at cr(Q). Then 
Y,jj is equivalent to the port Hamiltonian system Ex, as in \4>ty an d 14-$ , *- e -, Eg preserves the port Hamiltonian 
structure of the system {!]], if and only if H = B*Y*(YQ ¥*) . □ 

Proof. Proof of Statement (1). First we prove the necessity. Let n satisfy (|2|). Since —J = J* and R = R* , H* satisfies 

U*(-J- R) + L*B* = S*U*. (48) 

Postmultiplying (|4"8"|) by QH and adding it to © premultiplied by U*Q, yields 

R = --[(JTQI1S - GL) + aTQHS' - GL)*], (49) 

with R as in (j43|) . Postmultiplying (|48| by QH and subtracting it from ([2]) premultiplied by n*Q, yields 

J = i[(n*Qns' - gl) - (tl*qus - GL)*], (50) 

with J as in (J43J) . Hence (J-R)Q(H*QH) = (n*QIL)(S-GL), B = (n*Qn)G, with n*Qn invertible, since Q is assumed 

invertible, proving the claim. 

The sufficiency is a direct consequence of Lemma [A. 2 1 with P = (n*Qn) _1 and CH — B*QH. 

Proof of Statement (2). First we show the necessity. Let T satisfy ([3]). Since —J = J* and R — R* , T* satisfies 

T*Q* = Q{-J-R)T* + QBTZ*. (51) 

Postmultiplying (|5"Tj) by Q~ X T* and adding it to ([3]) premultiplied by TQ _1 , yields 

R = --[(QTQ _1 T* - KB*r*) + (QTQ- X T* - TZB*T*)*}, 

with R as in (|4"4"]) . Postmultiplying ([5T|) by Q^^^T* and subtracting it from ((3]) premultiplied by TQ _1 , yields 

J = -[(QTQ- l T* - TZB*T*) - (QTQ- l T* - TZB*T*)*}, 

with J as in (l44l). Hence 



(J - R)Q = Q- KH, B = TB, 

with Q = (TQ _1 T*) _1 invertible, since Q is assumed invertible, proving the claim. 

The sufficiency is a direct consequence of Lemma [A. 2 1 with P = TQ _1 T* and P~ l BT = H* . □ 

If S and L are fixed, then En is the unique model from the class Eg that matches the moments of (JTJ) at a(S) and 
preserves the port Hamiltonian structure. Assume now that S is fixed and let L = [li 1% ■■■ l v ), h € C, i = l,...,v, 
be such that the pair (L,S) is observable. Then the solution of the Sylvester equation $Z§ is given by a matrix n(L), 
yielding a class of reduced order port Hamiltonian models E n (L) defined by f4"2"]) with J(L), R{L), Q(L), B(L) as in (J4*3|) . 
Note that the input output behaviour is not affected by the choice of l\, ..., l v , i.e., all models parameterized in L have 
the same transfer function. However, since the port Hamiltonian structure is a state-space property, the parameters li, 
i = 1, ...,i/ can be used to enforce additional structure, e.g. the Hamiltonian defined by Q, or the dissipation R, have a 
desired form. Similarly, let H = [n ri ... r„]*, r* € C, i = 1, ..., u, be such that {Q,TZ) is controllable. Then the solution 
of the Sylvester equation ([3J is given by a matrix Y(1Z), yielding a family of reduced order port Hamiltonian models Ex(-r.) 

defined by equation (|42|) with J(1Z), R(7l), Q(1Z),B(1Z) as in (|44|) . All models parameterized in 1Z have the same transfer 
function and the parameters r^, i = 1, ..., v can be selected such that the models meet additional constraints. 



'Two minimal systems described by state-space equations arc called equivalent if they have the same transfer functions. 
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In the MIMO case, U £ C m , i — l,...,z/ can be chosen such that the right tangential interpolation conditions ([§]) are 

satisfied. Similarly, n G C lxp , i = 1, ..., v can be chosen such that the left tangential interpolation conditions (fTU)) are 

satisfied. 

Corollary 3. Consider a MIMO, port Hamiltonian system (QP, with the input u(t) £ W" 1 and the output y(t) £ M. m , i.e. 

B £ C nxm and C G C mx ™. Let L = [h ... Q € C mxi/ and ft = [r* ... r*]* G C xm . Then the following statements hold. 

1. The MIMO system En described by equations CM) and J^ffi ), ura'£/i input u(t) and output ip(t) £ M m , satisfies the 
right interpolation conditions |PJ). 



2. The system Ex described by equations CM) and C3], with input u(t) and output ^(i) G K m , satisfies the left 
interpolation conditions Iil0\) . D 

Theorem[5]offers a way to find a reduced order port Hamiltonian model, from a reduced order model that achieves matching 
of moments of the given (port Hamiltonian) system, by selecting the parameter G. Let ^ be a reduced order model and 
let P be such that S*P + PS < WQBL + L*B*QIi. Then, according to TheoremEU there exists G such that the model 
is passive, i.e., PG = H*QB. Selecting P = n*Qn yields the parameter G which identifies the port Hamiltonian reduced 
order model that achieves moment matching. Based on Lemma IA.21 similar arguments hold for the case of the system 
E T , with P = 1Q- X H*. 



Lemma IB . 1 1 implies that the result from Theorem [S] is equivalent to the result from Theorem [31 i.e., the class of reduced 
order models obtained by Krylov projection and the class of reduced order models obtained by time-domain moment 
matching are equivalent. In detail, let Eg and Ey be two reduced order models of fT}. Then selecting T = n*Qn 
yields U*QU(S - GL) = W*{J - R)QVU*Qn, n*QnG = W*B and B*QU = B*QVU*Q\l 1 which shows that one port 
Hamiltonian model can be obtained from the other via a coordinate transformation. Note that from an input-output 
point of view, both models exhibit the same behaviour. Similar arguments hold for the case of the system Ex, with 
T = YQ~ 1 Y*. Hence, from a computational point of view, Theorem [5] offers a way to compute a reduced order model En 
(or Ex) that achieves moment matching and preserves the port Hamiltonian structure. 

Algorithm 1. (Computation of a port Hamiltonian reduced order model that matches a prescribed number of moments 
of a given port Hamiltonian linear system, at a set of finite interpolation points.) 

1. Select S and L such that the pair (L,S) is observable. 

2. Use any efficient algorithm (e.g. Iterative Rational Krylov) to compute V depending on the eigenvalues of S. 



3. According to Lemma TB. 11 set 11 = V. 

4. Compute the class of reduced order models Eg as in ([6]). 

5. Let G = (Tl*QTl)~ 1 Tl*QB and compute the port Hamiltonian model E n as in (g2j and ([4"3"]l. □ 

The outcome of Algorithm [1] is a port Hamiltonian model that matches the moments of ((T|) at o~(S), parameterized in S 
and L. Hence, following arguments from Section |2~31 a suitable selection of the interpolation points, at step 1, such that 
the model approximates ([T]) in the H2 norm more accurately, is given by the matrix S such that <r(S) = {— Ai, ..., — A„}, 
where Ai G a ((J - R)Q). 

4.2 Markov parameter matching with preservation of port Hamiltonian structure 

In this section we solve another particular instance of [TJ for the case of matching a set of prescribed Markov parameters. 
Problem 3. Consider the linear, port Hamiltonian, SISO system ([I} and the observable pair (L, S), where L G C lxy and 
S £ C 1JXl ' such that Assumption Q] holds. Find F eC w and/or G £ C v , such that the following properties hold. 

(ml) The family of systems En parameterized in F and/or G and H, described by equations ((27)) . with II G {n, II}, 
where n is the unique solution of ([20| and n is the unique solution of ([2~T|) match the moments of the transfer 
function K{t) of ([1]) at cr(S), i.e., conditions (|28| or ((29)) are satisfied. 

(m2) There exists G and H such that the family of systems En, described by equations (|2"T)) . preserve the port Hamiltonian 
structure. 

(m3) The computation of the family of systems En described by equations ([2~T1) avoids the explicit solution of the Markov 
parameters. □ 

Based on the solution of subproblcm (ml) provided in this section, we provide solutions to subproblem (m2), suitable to 
the scope of this paper. We give some insight into subproblem (m3) provided by Proposition [5] which allows for efficient 
computation of the families of reduced order models. Note that a similar problem can be formulated in terms of the results 
of Section 13.11 hence it is omitted. 



In this section we focus on Problem [31 i.e., find a reduced order model that matches a prescribed number of Markov 
parameters and preserves the port Hamiltonian structure. Consider the port Hamiltonian system ([TJ) with the transfer 
function K(s) = B*Q(sI — (J — R)Q)~ 1 B. Suppose Assumption [TJ holds. The moments of K{r) = K(1/t) at r = r* are 
in a one-to-one relation with B*QHS, where n is the unique solution of the Sylvester equation 

( J - R)QUS + BL = n. (52) 
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In addition, let II be the unique solution of the Sylvester equation 

( J - R)QILS + BLS = EL (53) 

The moments of K(t) = K(1/t) at r = t* are in a one-to-one relation with B*QYl. The first v Markov parameters of 
|T]) are the moments of K(t) for r* = 0. Assume there exists an invertible matrix P such that a reduced order model 
described by equations ([77)1 exists and the relations (1251) are satisfied. Furthermore, assume there exists an invertible 
matrix P such that conditions (l29l) are satisfied. 



At the same time, we compute the class of port Hamiltonian reduced order models that achieve moment matching based 
on Proposition |U The moments of K(t) = K{1/t) at r = r* are in a one-to-one relation with QTB, where T is the 
solution of the Sylvester equation 

QT(J -R)Q + KB*Q = T. (54) 

The first v Markov parameters of ([T|) are the moments of K(t) for r* = 0. 

Proposition 8. Let K(s) be the transfer function of the system {!]]. Let (L, S) be an observable pair and let (Q,1Z) be a 

controllable pair. Consider the port Hamiltonian system 

i = {J-R)Q£ + Bu, 

(55) 

V. = B*Qt;, 

with £(£) € K y . The following statements hold. 

1. Let 

J = Tl*QJQU, R = TTQRQn, 

(56) 
Q = (U*QU)-\ B = U*QB, 

with II = II, or II = H, where II is the unique solution of equation \52\) and Pi is the unique solution of equation 
Ii53\) . Let Sn denote a port Hamiltonian model described by equations \55\) and i56}) . Lf X/i ^ 1, for any A G u(S) 
and any /i G a((J — R)Q), then the system En is a port Hamiltonian, reduced order model that matches the moments 
ofK{ T ) =K(l/s) ata(S). 

2. Let 

J = T*JT, R = T*RT, 

(57) 
Q = (T*Q" 1 T)- 1 , B = T*B, 



with T the unique solution of equation {5J$ - Let Ex denote a port Hamiltonian model described by equations A55\) 
and J57| ). If Xfi ^ 1, for any A G cr(Q) and any [i G o~((J — R)Q), then the system Ex is a port Hamiltonian, reduced 
order model that matches the moments of K(r) — K(l/s) at o~(Q). D 

Proof. Proof of Statement (1). Let II = II be the unique solution of equation (P0)l . Let P = Q^ 1 = H*QH, then 

(J - r)qps + bl-p = il*q(j- i?)Qn(n*Qn^ 1 n*Qn5 + itqbl - n*gn = n*q[(j - r)qils +bl - n]. By duji 

we have that (J - R)QUS + BL-U = 0, hence (J - R)QPS + BL = P. Furthermore B*QP = B*QU(U*QU)~ 1 U*QU = 
B*QU. Then, system (|55| satisfies conditions (l28l) and the claim is proven. The proof is similar in the case II = EL 
The proof of the second statement follows similar arguments as the proof of Proposition \7\ hence it is omitted. □ 

Proposition |5] shows that there exists a port Hamiltonian model that matches the moments of a given port Hamiltonian 
system by direct computation of the matrix II. This computation is avoided by showing that the model En is a member 
of the class of reduced order models ([27| that match the moments of K(t) = K(1/t) at <j(S), for a particular instance 
of the parameter G and of the output H . Similarly, we also show that the model Ex is a member of the class of reduced 
order models (|57| that match the moments of K{t) = K(\/t) at a(Q), for a particular selection of H. 
Theorem 7. The following statements hold. 

1. Let system En, as in {27}) , be a reduced order model of system (QJ) that matches the moments of K(t) — K(l/s) at 
o~(S), where K(s) is the transfer function of (QJ). Then En is equivalent to a port Hamiltonian system En, described 
by equations i55\) and A56\) . if and only if G = Tl*QB and H — G*(II*Qn) _1 . 

2. Let system |ff7| j be a model that matches the moments of K(t) = K(l/ s) at a{Q), where K(s) is the transfer function 
of (Qp. If E = I, then \3ty is equivalent to a port Hamiltonian system Ex described by equations \53i) and J57| ), if 
and only if H = B*T*Q. □ 

Proof. Proof of Statement (1). The necessity follows from the fact that ([55)1 matches the moments of (HJ. Let II = n 
satisfy (|52|) . Postmultiplying the transpose of ((52|) by QH and adding/subtracting it to/from (|52l) premultiplied by U*Q, 
yields 

(J-R)S + S*(J + R) = L*B-BL, (58a) 

(J- R)S - S*(J + R) = 2n*gn - L*B - BL, (58b) 
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with R and J as in ([56")) . Furthermore, adding equations (|58ap and (I58b|) yields 

(j - p)Qn*Qns = n*Qn - pp 

with ITQII invertible, since Q is assumed invertible. Moreover, P*Qn(n*Qn)- 1 n*QnS' = B*QTIS J Hence (J55D is a 
model described by equations J5JJ with F = (J - R)Q, G = B, H = G*P" 1 and P satisfying P = Q" 1 = irQII from 
the matching conditions (l28l) . 

The sufficiency is proven in two parts. First we assume that S has full rank. Let (p?T|) be a reduced order model of (UJ), 
with G = WQB and H = G*(U*QU)~ 1 . System (J2ZJ) matches the moments of K(r), since there exists P = ITQn such 
that FIl*QIlS + Tl*QBL = n*Qn, Pn*QnS' = B*QUS. In addition, we have FTTQUS = n*Q{U-BL), which by flgg) 
becomes [P — ( J — P^QjrFQIIS' = 0. Since S 1 is assumed to have full rank, F = (J — R)Q, i.e. (P2T1) is a port Hamiltonian 

" Si " 



system described by equations (|5"5j) . If 5 does not have full rank, assume, without loss of generality that S 



S 2 



with S^ square, of appropriate dimensions and full rank, and L = [1 ... 0]. Let II = [0 II], with II (having appropriate 
dimensions) satisfying the Sylvester equation (J — P)QII52 + BS\ = II. Assume (p?T)) is a reduced order model with 
G = H*QB and H = G*(H*Q H) _1 . Since S2 has full rank, the arguments from the previous case are followed. Hence, 
applying the matching conditions ([31]) yields [P — (J — P^QJIFQIIS^ = 0, which leads to the claim. The proof is similar 
for the case II = II. 

The proof of the second statement follows similar arguments as the proof of Theorem [6l for P = I = Q~ l Q, where P is 
an invertible matrix uniquely satisfying the relation (|38|l . □ 



If r* = 0, the models E£j match the first v Markov parameters of ([T]) and preserve the port Hamiltonian structure of the 
given system. Furthermore, a direct application of Proposition [5] yields that Ey = Sri, with £y as in (jlip and (fT3 



Based on Theorem [7] and on Proposition [5] we propose the following procedure. 

Algorithm 2. (Computation of a port Hamiltonian reduced order model that matches a prescribed number of Markov 

parameters of a given port Hamiltonian linear system.) 

1. Select a non-derogatory matrix S, such that a(S) = {0, ..., 0}, and L such that the pair (L, S) is observable. 



Use any efficient algorithm to compute V as in (I40al) . 

Set n = V. 

Compute the class of reduced order models En as in (071) . 

Let G = H*QB and H = G*(H*QU)~ 1 and compute the port Hamiltonian model En as in (J55J) and ([56 



□ 



The outcome of Algorithm [2] is a port Hamiltonian model that matches the Markov parameters of ([T]) , parameterized in 
L. The parameter L cam be used to enforce further additional structure on the state-space realization of the reduced 
order model, increasing the physical meaning of the approximant. Note that if the result of Algorithm |2] is Ejj with fl the 
unique solution of (EJJ, then we obtain a i^-th order model that matches v + 1 Markov parameters, see also Example [5] 

Note that the results of both Theorem [6] and Theorem [Jj apply to the port Hamiltonian system in any coordinate sys- 
tem, resulting in equivalent state-space port Hamiltonian reduced order models, without the need to compute additional 
transformations. 



5 Illustrative examples 



5.1 SISO ladder network 



u = I 



y = Vcr 



Pi £i,0i P2 l 2 , 4>2 



Ci,qi 



C 2 ,q2 



R-3 



Figure 5: Fourth order ladder network. 

Consider the ladder network in Fig. [5J with C\, C 2l L±, L 2 , Pi, R 2 the capacitances, inductances, and resistances of 
the corresponding capacitors, inductors, and resistors, respectively. The port Hamiltonian representation of this system 
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is given by equations of the form ((T|), with x = [q\ <j)\ q 2 ^2]* an d 



J = 



0-100 
10-10 
10-1 
10 



R = di ag {0,R 1 ,0,R 2 + R 3 }, Q = diag<j -L,-L,-L,J- J>, B = [10 0]*. C<)i 

L-i ivi O2 Li 2 



Note that </> denotes the flux through the inductor L and q denotes the charge at the capacitor C . Assume C\ — C 2 7^ 0, 
L\ = £2 7^ 0, R\ = R 2 = -R3 7^ 0. The transfer function of the port Hamiltonian system ((59|) is 



#(*) 



L 2 Cis 3 + 3Lii?iCis 2 + (2Li + 2i?fd) s + 3i?i 
C^L 2 ^ + ^l^r^3 + (3(7^! + 2C\R\) s 2 + 5R^h~s + 1 ' 



Matching at finite interpolation points. The first two moments of (|59[) at are 770 = 3i?i and 771 = 2L\ — l'ZR\C\. 

" 1 " 



Let L = [1 0] and 5 







Note that 



n = 



3i?iCi Ci(£i - 13Cii? 2 ) 
.Li — 3-fti.LiGi 

2R l C l C X {L X - 10i? 2 Ci) 
-Li — oR\L\G\ 



A reduced order port Hamiltonian model that matches the moments 770 and 771 is given by (|42p . with 



J = 







2Ii 
-2ii 



R — R\ 
1 



3 

-13.RiCi 



10Z,f - i? 2 Ci(HL 2 - 44i? 2 ii - 16i?f) 



-13i?iCi 
59RlCf 

5L 2 -i? 2 Ci(38ii-269i?f) 
59i?fC*i 



59i??Ci 



B 



3i?i 
2Ii - 13i? 2 Ci 



(60) 



The transfer function of the reduced order model is Ku(s) 



6 s 2 + f s + § 



, with a, b, c given by 



a = Rfd(mCl + 28i? 2 LiCi - 7L\) + 8Lf, 
& = Ci(10L? - i?fCi(lli 2 - 44i? 2 Li - 16Rf)), 
c = J RiCi(40i?iC* 2 + 15L 2 + 4i? 2 iiCi), 
d= \2R 1 {L\ + 2R\Cl), 

e = 4(Ll + 2RtC 2 ). 

Comparing the port Hamiltonian model En in (|60[) with the model Ey in (jlll) yields En = E_y, i.e., n = — V, where 
V spans a Krylov subspace of the given port Hamiltonian system given by (l59|) (see also Proposition [2]) . Note that this 
yields a reduced order model with the transfer function Kn(s). 



Let d = 1,C 2 = 2,L X = L 2 = 1 and i?i = R 2 

observable if and only if /1 7^ 0. Note that 

n(h,i 2 ) = 



R 3 = 1. Furthermore, let L = [h l 2 ], h el, l 2 G K. The pair (£, S) is 



3/i 



h 



h 



h 



3^2 — h l 2 — 3/l ?2 ~ 2^ ^ 2 _ ^1 

The class of port Hamiltonian models, parameterized in /1 and l 2 is given by 



J(h,h) = 



2Z? 
-2/ 2 



, R(h,l 2 ) = 



3/ 2 



3/ 2 /i - lllj 



3/ 2 /i -lllj Zll - 22l x l 2 + All\ 



Q{hM) = ^l 



1 
3l7f 



26/1 - 164ZiZ 2 + 261/ 2 2/i(41/i - 13/ 2 ) 



2/i(41/i - 13/ 2 ) 



26/f 



B(h,l 2 ) = [3/i 3/2-9/1]* 



For l 2 = ji/i, we obtain the subset of reduced order models with the following properties: they match the first two 
moments of (|59p at zero, preserve the port Hamiltonian structure of the model and have diagonalized Hamiltonians. 
For Z2 = ^r/ii we obtain a subset of port Hamiltonian reduced order models with diagonal dissipation matrix. All the 

parameterized models have the same input-output behaviour described by the transfer function Kjj(s) — 31f , 2 , 4 5s li 2 ■ 
Note that the physical meanings of the states of the approximant are more difficult to recover using only the parameters 
h, i = 1, 2. However, the selection of the interpolation points Sj, i = 1,2 (as free parameters), can bring additional insight 



into the physics of the approximant. 
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Q{h,k,h) 



' -If 

l\ 



h(h-k) 

-Ml + Zi Z3 — l 2 + 'l'3 





R(h,k,h) 





l\ h(h-h) 

ohih-h) {W-hf 



1 

2lf 



gii(Ji,J2,fa) qi2{h,h,h) qi3(h,k,h) 
qi2(h,h,h) q22(h,h,h) q2a(h,k,h) 
qi3(h,h,h) qs3(h,k,h) qsa(h,h,h) 



, B(i l ,k,h) = [hkh-hY 



(61) 



"Dual" family of port Hamiltonian reduced order models. Let 1Z = [r\ r 2 ]* , r\ € R, r 2 G K.. The pair (Q, 72.) is 
controllable if and only if n ^ 0. Solving (|3]l. we obtain 



T(ri,r 2 ) 



3ri — 7"i 2ri — ri 

3r 2 - ri 3ri — r 2 2r 2 — 7r\ 4ri - r 2 



The class of port Hamiltonian models, all with transfer function Kt{s) = 302J27+6 ' P arame terized in ri and r 2 is given 

by 



J(ri,r 2 ) ■■ 

Q(n,r 2 ) 



-2rf 

2rf 

1 



, R(r!,r 2 ) 



3r 2 3r 2 ri - llrf 

3r 2 ri - llr? 3r| - 22rir 2 + 41rjf 



32rf 



204rf - 124rir 2 + 19r| r x (62ri - 19r 2 ) 
rx(62ri - 19r 2 ) 



19r? 



B(ri,r 2 ) = [3n 3r 2 - 9n]* 



For r 2 = firi we obtain a subclass of port Hamiltonian reduced order models with diagonal Hamiltonians. For r 2 = -g-ri 
we obtain a subclass of port Hamiltonian reduced order models with diagonal dissipation matrix. 



Matching at infinity. Let L = [h k h]*, h € R, Z 2 € R, Z 3 € R. Note that 

n(j 1 ,i 2 ,fe) = 



Zl 


k 


h 


-ill 





h 


h 


-Jl 










Zl 














The family of port Hamiltonian models, parameterized by Zi, Z 2 , Z3, is given by equation (|6ip with 

«n(*i> *2, fe) = 5/fZl - 2/^3^ +4- 2llh + il\ - 21% 

+ Z 3 Zj^ — 2Z 1 Z 2 + 2I2I3I1, 
qi2(h,h, h) =lf~ 4Z?Z 2 + ZfZaZs " 1% + 2l\l\ - l\k, 
qi 3 (h,k,h) = lUl2-hk + ll-hh) ) 

q22{hMM)=l\{ll~2l 1 k + ? > ll), 

q2z{h,k,h) = l\{h-hh), 
qza{h,k,h) = If 

The input-output behaviour of the family of models described by (|61[) is given by the transfer function Ku(s) 



s^+s+2 



Let now L 2 = [k 0], Si = [1 0], 5 2 
particular case is 



1 




L = [h L 2 ] and 5 



Si 
S 2 



(s 2 +s+3) • 

Consider now equation (|30[) . which in this 



n = 0, (J - i?)Qns 2 + BiiSi + bl 2 ^2 = n. 

The second and the third Markov parameters of (|59p are in one-to-one relation with B*QH, with n(Zi, Z 2 ) = 



(62) 



(1 000 
k h 

A family of reduced order port Hamiltonian models, of dimension two, that match the first three Markov parameters is 
described by equations ([55]) . with 



J(h,k) = 



-l\ 
l\ 



, R{h,k) 





oil 



Q(Ji.fc) = 7J 



*? + i| 
-w 2 



-Z1Z2 

I 2 



B(h,k) = [hk}*. 



s+l 



S 2 +S+ 



j . Note that selecting 



The input-output behaviour of this family of models is given by the transfer function Ku{s) 

l\ = 1 and k = 0, yields Q(Zi) = 7 and Z?(Zi) = [1 0]*, rendering the reduced order model a ladder network as in Fig. [6] 

with the parameters L = C = (7?i + 7? 2 ) = 1. Furthermore, let 



Q = 





1 



, K = 



r 2 
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u = I 

— > — 



j, = V s : 



Ri 



L, 



:C,q 



Ri 



Figure 6: Second order ladder network. 

with ri ^ 0, r 2 € K. Let T be the unique solution of ©. The first two Markov parameters of K(s) are in one-to-one 
relationship with Q(QT + TZB*Q)B = [0 n]*. Consider J, R and Q as in (57]), i.e. 



J = T JT* 

where 



3r? 
-3r? 



R = TRT* 



rf ri(r 2 -ri) 

n(r 2 -n) (n -r 2 ) 2 



Q=(TQ- 1 T*)- 1 = - T 
6r l 



4rj-2r 1 r 2 +r 2 : 


ri— r 2 


ri-r 2 
ri 


1 



T = 



o -n oo 

-n ri — r2 2ri 



is the unique solution of ((2). Furthermore, solving (|31?]) yields 

P = diag{V3, V3}, B = 







Then a port Hamiltonian reduced order model that matches the first two Markov parameters of the given system Ex, 

8+1 

s 2 + s+3' 



described by equations (|55|) and (|57|) with the transfer function Ky(s) — s+1 



5.2 MIMO single machine infinite bus system 



synchronous 
machine 



rrn 



transmission line 



infinite bus 



Figure 7: A single machine connected to an infinite bus through a transmission line. 

Consider the model of a single machine connected to an infinite bus (SMIB) useful in the analysis of power systems 
stability, where a power system is modelled as the interconnection of a large number of such systems. The machine under 
consideration has one field winding, three stator windings, two g-axis amortisseur circuits and one d-axis amortisseur 
circuit, all magnetically coupled consisting of electrical and mechanical equations, see e.g. [32lH]. 



The nonlinear port Hamiltonian model of the SMIB is given by equations of the form 

x = (J(x) — R)Qx + Bu, y = B*Qx, 

with 



J(x) 



R — diag{i? a -I- Re, R a + Re, Rfd, Rid, Ri q , R 2q , Kd\, 

and Q = diag < L _1 , i \, where 

• 5 and u are the angle and the angular velocity of the rotor, respectively, 

• Re and Xe are transmission line resistance and reactance, respectively, 



(63) 






loX e 





*<r 




sin S 


-lqX e 








-*d 




cos<5 
















1 














, B = 







































-*« 


*rf 










1 



22 



• ^d and $ q are the stator fluxes in the ci-axis and g-axis, respectively, 

• R a is the stator resistance, Rfd is the field circuit resistance, Ri, i £ {Id, lq, 2d} are the amortisseurs resistances, 

• Kd is the damping constant and j is the inertia of the rotor, 

• L > is the inductance matrix, see, e.g., [39] for further details. 

The state is x — [** w]*6t 7 , where 

with 

• ^Jds = ^d + XE%d and ^s = ^ + Xsiq, where id and i q are the stator currents, 

• */ d the field flux, 

• ^i, is {Id, lq, 2d} the rotor fluxes due to the amortisseurs. 

The input is u — [Et e/d T m ] £ K 3 , where E}, is the infinite bus voltage (interconnection variable), ejd is the field voltage 



(control variable) and T m is mechanical input power. The (passive) output is y = B*Qx = [lb ifd w m J 
infinite bus current, ifd is the field current and co m is mechanical output angular velocity. 



, where h is the 



Linearising the model (1631) around an equilibrium point x* = [$* p*]* yields a system described by equations of the form 
([1}, with J = 0, defined by (-RQ, B, B*Q), which is minimal and passive. Note that, we solve the tangential interpolation 
conditions ©■ To this end, we apply Algorithm Q] to the study case defined by the equations (j6"3"|). We have the following 
parameters for the machine, taken from e.g. [32] 



0.22 0.01 0.01 

0.219 0.009 0.009 

0.01 1.825 1.660 

0.01 1.660 1.8313 

0.009 0.009 

0.009 0.009 0.134 



B = 



1 

10 

0.7071 0.7071 



R = diag {0.031,0.031,0.0006,0.0284,0.00619,0.023638,10} and j = 6. We consider an equilibrium point which is 
asymptotically stable. Linearising around this equilibrium point we obtain a minimal asymptotically stable linear realiza- 
tion with S as a parameter, see, e.g., [25]. Let S = diag {0.055, 0.01, 1.667, 0.0021}, h = [1 0]*, l 2 = [0 1 0]*, l 3 = [0 1]*, 
I4 = [1 1]* and L = [li I2 I3 I4]. Note that (L,S) is observable. Using any numerically efficient algorithm compute 
U = V = [{si + RQy 1 Bh (s 2 + RQ)~ 1 Bl 2 (s 3 + RQ)' 1 Bl 3 {s A + RQ)~ l Bl A ]. The fourth order linear system that matches 
the moments of (—RQ,B,B*Q) at (L,S), in the sense of satisfying the right tangential interpolation conditions ^ is 
(S-GL,G,B*QU), with 



S-GL 



-1.6667 -0.0048 
-0.0000 -0.0081 
0.0000 -0.4825 
-0.0000 0.0047 



-0.0004 -1.7220 
-0.0002 -0.0002 
-0.0875 -1.7545 
0.0004 0.0025 



, G 



1.7217 0.0048 0.0004 
0.0000 0.0181 0.0002 
-0.0000 0.4825 1.7545 
0.0000 -0.0047 -0.0004 



Note that II is the unique solution of @. Using Theorem [51 the port Hamiltonian state-space representation of (S 
GL, G, B*QU) is given by a system (g2J| with 



R = P(S - GL) = 



0.0937 -0.0000 

0.0000 71.6093 

-0.0000 -0.2859 

0.0967 -65.0309 



0.0000 0.0967 

-0.2859 -65.0309 
0.0525 1.7312 

1.7312 135.8342 



B = PG = 



0.0968 -0.0000 0.0000 

0.0000 162.9925 -73.4077 
-0.0000 -0.7246 1.9466 

0.0999 -104.9205 172.7502 



with P = (rTQlLr 1 = Q- 1 , where 



Q 



17.7849 -0.0000 
-0.0000 0.0001 
0.0111 0.0051 
-0.0001 0.0000 



0.0111 -0.0001 
0.0051 0.0000 
2.0518 -0.0108 
-0.0108 0.0001 



Note that (S — GL,G, B*QH) is a passive system. For an accurate approximation, following the arguments of Remark 
[TJ the interpolation points chosen were approximations of the mirror images of the reduced order model. The initializing 
choice were the poles of the fourth order balanced truncation, known to have a good approximation error. 

Figure [S] shows the Bode plots of the elements of the transfer matrix that correspond to the transfer between the passive 
inputs and outputs respectively, i.e., the transfer function from the bus voltage Ef, to the bus current %, as in Fig. 8(a) 
the transfer function from the field voltage Efd to the field current Ifd, as in Fig. |8(b)| and the transfer function from the 
torque T m to the angular velocity w m , as in Fig. 8(c) respectively. The plots show both the original and the approximated 



responses, depicted by solid line and starred line, respectively. 
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Frequency (rad/s) 

(c) 

Figure 8: Bode plots of the transfer function between the bus voltage E\, and the bus current if, (a), from the field voltage 
Efd to the field current Ifd (b) and from the torque T m to the angular velocity uj m (c). The solid line represents the 
evolutions of the 7th order SMIB and the starred line represents the evolution of the reduced 4th order model 
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6 Conclusions 

In this paper we have discussed the solutions of the port Hamiltonian structure preserving model reduction Problem [31 
based on time-domain moment matching at a set of both finite and infinite points. In the case of matching at finite 
interpolation points we have first obtained the reduced order port Hamiltonian model that matches the moments of a 
given port Hamiltonian system (Proposition [7]). Furthermore, we have characterized all the reduced order models which 
match the moments of the given port Hamiltonian system and preserve the port Hamiltonian structure (Theorem [S]). We 
have obtained families of state-space parameterized, reduced order port Hamiltonian models that approximate the given 
port Hamiltonian system, all models having the same transfer function. The state-space parameters allow to enforce 
additional constraints on the structure and/or state-space realization and, in the MIMO case, they have been used to 
solve the tangential interpolation problem. We have given a possible procedure to compute the family of port Hamiltonian 
approximations (Algorithm [IJ . 

We have also studied the problem of Markov parameters matching, extending the time-domain moment matching results to 
the case of interpolation points at infinity. We have defined the moments of a class of linear, descriptor systems, associated 
to a given linear system, in terms of the unique solutions of Sylvester equations and their dual counterparts (Propositions^ 
and[¥]). In particular, the Markov parameters of a given system are the moments of a descriptor realization associated to the 
given transfer function at zero. Furthermore, we have related the moments to the well-defined steady-state response of the 
descriptor realization driven by/driving signal generators (Theorems 2] and [5]). We have obtained families of parameterized, 
descriptor reduced order models that match a set of prescribed moments of the descriptor realization associated to a given 
linear system (Proposition [5]). In particular, matching at zero has yielded classes of reduced order models that match 
the Markov parameters of the given linear system. Finally, applying these results to linear port Hamiltonian systems, we 
have solved Problem [31 yielding families of state-space parameterized, reduced order port Hamiltonian models that match 
the Markov parameters of the given port Hamiltonian system (Proposition [51 Theorem [7] and Algorithm [5]) . Finally, the 
examples proposed in Section [5] have illustrated the aforementioned results. 

For future work, the nonlinear extension of the results in this paper is a goal. Furthermore, it is very important to 
determine how many variables are retained in reduced order model, the relation of the variables between two systems, and 
the physical meaning of each variable and we will address this issue in the future. 
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A Preliminaries for the proof of Theorem [7] 

Theorem A.l. JfJjj The family of v-th order models Eg as in (0|) contains a passive system if and only if there exists a 
symmetric and positive definite matrix P such that 

S*P + PS<U*QBL + L*B*QIL, (A.l) 

where n is the unique solution of (0). □ 

Lemma A.l. A family of models (HI) contains a passive model if and only if there exists P = P* > G R"*" such that 
PQ* + QP<KB*T* -TBTZ*. D 

Proof. The proof follows immediately from applying the Kalman-Yakubovitch-Popov (KYP) lemma to system ([7]) (see, 
e.g., [33] for an original version and 031 051 S3] for extended versions). □ 

The next result shows how to obtain a port Hamiltonian system from the families of models Eg and £# , respectively, 
described by equations © and (J7J, respectively. 
Lemma A. 2. The following statements hold. 

1. Let Tiq, as in ftjjty, be a passive reduced order model of the system (QJ) and let P satisfy L4. 1\) . Then the matrices 
J = |[(5- P- 1 U*QBL)P- 1 - P~ 1 (S - P- 1 \TQBL)*\, R = -\[{S - p- 1 U*QBL)p- 1 + p-\S - P^WQBL)*}, 
Q = P and G = P~ 1 H*QB are such that Eg is a port Hamiltonian model En described by equations of the form 
g§ and g3p. 

2. Let Ejj, as in II, be a passive reduced order model of the system (Q]) and let P be as in Lemma \A.1\ Then the 
matrices J '= \[P{Q-1IH) - (Q-KH)*P], R = -\[P{Q -TZH) + (Q -HH)*P], Q = P' 1 and H = (P^TB)* 
are such that T,h is a port Hamiltonian model Ex described by equations of the form \4'2fy and \4$ . □ 

Proof. The first statement is identical to [301 Theorem 3] . The proof of the second statement follows directly from the 
application of Lemma IA.1I □ 

B Relations between the Krylov projections and the solution of the Sylvester 
equations 

Lemma B.l. f7/ The following statements hold. 

1. Consider the matrix H, solution of the Sylvester equation (0) and the projector V defined by equation H3\). There 
exists a square, non-singular, matrix T £ C l/Xv such that II = VT. 

2. Consider the matrix T, solution of the Sylvester equation {3|) and the projector W defined as in Theorem^ There 
exists a square, non-singular, matrix T G fry*-" suc h that T = TW . □ 
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